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Abstract 

We present a theory of the pair distribution function g{z) and many-body effective electron- 
electron interaction for one dimensional (ID) electron liquid. Our approach involves the solu- 



tion of a zero-energy scattering Schrodinger equation for W g(z) where we implemented the Fermi 
hypernetted-chain approximation including the elementary diagrams corrections. We present nu- 
merical results for g{z) and the static structure factor S(k) and obtain good agreement with data 
from diffusion Monte Carlo studies of the ID system. We calculate the correlation energy and 
charge excitation spectrum over an extensive range of electron density. Furthermore, we obtain 
the static correlations in good qualitative agreement with those calculated for the Luttinger liquid 
model with long-range interactions. 
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1. INTRODUCTION 



Electron-electron interaction is known to produce the most pronounced effects in one 
dimensional electron liquid (ID EL) systems, where a strongly correlated state appears [lj]. 
The effects of interactions in higher than one dimensional systems has been masterfully 
explained by Landau's Fermi liquid theory [2]] dealing with low-lying excitations in a system. 
Landau called such single-particle excitations quasiparticles and postulated a one-to-one 
correspondence between the quasiparticles and the excited states of a non-interacting Fermi 
gas. However, for the ID EL any individual excitation becomes a collective one. This 
collectivization of excitations is obviously a major difference between the ID and higher 
dimensions electron liquids. It clearly invalidates any possibility to have a Fermi liquid 
theory in ID EL. One of the most interesting characteristics of the system is the spin-charge 
separation seen in a series of remarkable experiments j^]. Furthermore, according to the 
theorem proved by Lieb and Mattis jj], theground state of a ID EL is paramagnetic. The 
recent quantum Monte Carlo calculations [5] consonant with the Lieb and Mattis theorem 
predicts no Bloch instability [6] for the ID EL. 

Many of the electron-electron interaction effects become increasingly important as car- 
rier density and dimensionality are reduced and the homogenous electron liquid provides a 
primitive model for their study. The crucial role in the theory is played by the particle pair 
distribution function g(z). In our previous works we presented a theory of the pair distri- 
bution function and other ground-state properties of the 3D and 2D electron liquid 0, S] • 
The theory was based on a Fermi hypernetted-chain approximation (FHNC) which rep- 
resents a direct generalization to the Fermi systems of the well known hypernetted-chain 
Euler-Lagrange (HNC/EL) approximation. One of the important feature of theory based on 
HNC/EL or FHNC is that the theory sums not only all ring and ladder diagrams exactly, 
but also mixed diagrams in a local approximation [9]. It means the FHNC theory is not a 
Fermi liquid theory and therefore we able to apply it for studying many-body effects on ID 
EL. As we shall see, much higher sophistication is needed to attain a quantitatively useful 
theory in ID EL. We shall have to dwell on terms beyond the FHNC/0, which includes 
elementary diagrams and the three-body Jastrow-Feenberg correlations. These effects have 
been studied theoretically in boson fluids for the ID system using the HNC/EL formalism 
by Krotscheck et al [lO ]. 

The main purpose of the present work is the formulation of a practicable theory for charge 
correlations in the paramagnetic ID EL and comparison of its numerical predictions with 
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the available diffusion Monte Carlo data 



ll| . The theory is based on variational FHNC 
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method made by simplifying the standard FHNC method to have only one Euler-Lagrange 
Schrodinger equation. The computational time typically needed to obtain g(z) at each value 
of density is a few minutes on a PC with a Pentium IV/3.4 GHz processor. There are many 
theoretical works addressing the ground-state properties of the quasi- and strictly [13j, \14 \ 
ID EL but an accurate charge correlations is still missing. 

In this work, we numerically study the many-body effects for ID EL by calculating the 
pair distribution function and electron-electron effective interactions. 



THEORY 



We consider a ID EL as a model for a system of electronic carriers with band mass 
m in a semiconductor heterostructure with dielectric constant e. The resulting effective 
ID potential is readily shown 15| to be v(z) = (e 2 /e)(v / 2/26) exp[z 2 /4b 2 ] erfc[\z\ /2b] 
with Fourier transform v(k) = (e 2 /e) n exp[b 2 k 2 } E\\b 2 k 2 ] where the exponential integral 
function, Ei(x) is defined as J x °° e~ u du/u. The Fourier transform is defined according to the 
general expression FT[F(r)] = n J drF(r) exp (ik ■ r). Here n is the total average density. 
The above form of the bare potential exhibits the typical ID behavior, v (k) « ln(kb) as 
k — > 0, and the 3D behavior, v(k) « l/(kb) 2 as k — > oo. At zero temperature there are 
only two relevant parameters for a disorder- free, homogeneous ID EL in the absence of 
quantizing magnetic fields, namely the quantum wire width b and the usual Wigner-Seitz 
density parameter r s = (2na* B )~ 1 , where a* B = h 2 e/(me 2 ) is the Bohr radius in the medium 
of interest. We assume that only one subband is occupied and neglect any contribution from 
higher subbands. This approximation leads to r s > irb/A. 

With the zero of energy taken at the chemical potential, the formally exact Euler- 
Lagrange equation for the spin summed pair distribution function g(z) reads 16l. Il7| 



-— + v(z) + v B (z) + v F {z) 
m az z 







(1) 



Here, v(z) is the ID EL potential and the Bose-like potential Vb(z) contains the effects of 
correlations and by itself would determine g(z) in a Bose fluid. The details of the Bose-like 
potential have been discussed by many authors, for instance see Ref. [181] . We can write it 

as 

V B (z) = W ind (z) + Av e l e (z) , (2) 
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where the "induced interaction" is ll 

h 2 k 2 

Wi n(l {k) = — — [S(k)-l}-v ph (k). (3) 

A momentum space formulation of the Euler equations equivalent to Eq. (jrj can be written 
in terms of the structure factor S(k), 

S(k) = - l L= , (4) 
1 + fii^2 v ph(k) 

and the FHNC-Euler-Lagrange theory supplements Eq. (TjO) with a microscopic theory of the 



particle- hole interaction 18] 



Vph(z) = g(z)[v(z) + Av e i e (z) + v F (z)] H 



h 2 ' 



+ [g(z) - l)w ind (z) . (5) 



m 

From the theoretical point of view it is important to remark that in FHNC-type calculations 
at strong coupling the Bose-like vb{q) interactions should be corrected by the addition of 



three-body correlations and elementary-diagrams (or "bridge functions") contributions 16, 
18 . I9I . 2o| . The Av e i e (z) is a term arising from elementary diagrams and triplet correlations. 
We only consider the fourth-order elementary diagrams and triplet correlations, Av e i e (k) = 
Wg i (k) + vjg (k), which hereafter, we call it as FHNC/4+triplet acronym. Note that the 
approach is reduced to the Fermi hypernetted-chain approximation, FHNC/0 when the 
corrections of the Bose-like interaction are omitted, Av e i e (k) = 0. The contribution from 
the low-order elementary diagrams to the effective Bose potential is 

where e^q) is given by a twofold integral (in ID) over momentum space. 

£ 4(9) = 2^ J ^[S(p)-l][S(q')-l][S(p + q')-l][S(p + q )-l] [S(p + q +q) - 1] , (7) 

and the contribution of three-body correlations is given by the integral, 

w { B } {k) = J dqS(p)S(q)u 3 (q,p, k) {v(q,p, k) + [e(p) + e{q)\ u 3 (q,p, k)} (8) 

where p = —(q + k), e(k) = h 2 k 2 /[2mS(k)] and, with the definition X(k) = 1 — S' _1 (/c), we 
have v(q,p, k) = (H 2 /m) [kpX{p) + kqX(q) + pqX(p)X(q)] and 

_ (h 2 /2m) [kpX(p)X(k)+pqX(p)X(q) + kqX(q)X(k)] 
U3[q ' P ' ] ~ e(k) + e(p) + e(q) ' [J) 



With regard to the Fermi term vp(z) in Eq. ([I]), we adopt the same criteria in determining 
its form as in Ref. An important requirement is that Eq. ([T]) should give the exact 

fermion distribution function when one goes to the weak-coupling limit r s — > 0, where in 
this limit g(z) becomes the Hartree-Fock pair distribution function, g™(z). The Fermi term 
in the scattering potential is then determined by the Hartree-Fock structure factor S HF (k), 

2 



h 2 

v F (k) = —FT 
m 



d 2 y/g HF (z)/dz 2 



^ h 2 k 2 \s HF (k) - i 



4m 



S» F (k) 



[2S HF (k) + l}-Av ele (k)\ 

(10) 

The Hartree-Fock static structure function is given by S™(k) = (k/kp)0{'2kF—k)+8(k—2kF) 
where 6 is the Heaviside step function, kp = 7r/Ar s a* B is the Fermi wave vector and g HF (z) 
is given by 

*"'M--l£f Pi) 



The rational behind Eq. (fTUj) . is as in Ref. 17| . (i) the first term on the RHS ensures that 



the Hartree-Fock limit is correctly embodied into the theory, (ii) the second and third terms 
ensure that the Bose-like scattering potential is suppressed for electrons at weak coupling [tJ . 
We introduce at this point the relation between g(z) and S(k): 

g(z) = l + ^ J cos{kz)[S{k)-l]dk . (12) 
2.1. Sum rules and limiting behaviors 

In this section we show that the pair distribution function obtained from the theory 
presented in Sec. II satisfy three exact properties. These are the behavior of the static 
structure factor at small momentum, the charge neutrality and the Kimball's cusp condition. 

The asymptotic behavior of the effective potential v p h(z) is first obtained from Eq. (jSJ). A 
careful analysis of this equation shows that v p h{z) — > v(z) + Vf(z) for z oo, and therefore 
in Fourier transformation we have v p h(k) v(k) + h 2 k 2 F /m for k — > 0, with v(k) as the 
Fourier transform of bare potentail the ID EL. The corresponding asymptotic behavior of 
the structure factor is obtained from Eq. (j4j), which we have 

hk 

S{k) - _^t= , (13) 
y / 4mv(k) 

for k — >• 0. The analytical behavior of correlations for the ID EL with long range potential 
interaction has been obtained by Schulz jgji using the bosonization techniques applied to 
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an effective ID hamiltonian with linearized kinetic energy expression. According to that 
model, S(k) behaves like oc kj *J\ln(k)\. It is clear that our scheme gives 

lim k ^S{k) = ■ (14) 

This result agrees with the bosonization findings in ID. The form of Eq. (I14p immediately 
ensures that the charge neutrality condition is satisfied. The charge neutrality condition, 

n j [g{z) - \}dz = -1 , (15) 

is equivalent to S(k — > 0) = 0, a result obtained from Eq. (fl4l) . We have numerically 
evidence that the charge neutrality condition is fulfiled for our total effective potential in 
Eq. (CQ). On the other hand, in the small— r limit the bare repulsion potential dominates 
over the induced plus Fermi potentials. Therefore, it follows from the differential Eq. (CQ) 
that 

0(* 0) = 0(0) (l + ^ *) • (16) 

This is a similar result of the Kimball's cusp condition |22| for 3D electron gas. Linear 
small— z behavior together with general properties of Fourier transformations implies that 
at large k the static structure function S(k) — 1 approaches zero asymptotically in proportion 
to -1/k 4 . 



2.2. Charge excitations in long-wavelength limit 



From the knowledge of static structure factor we can find the behavior of the charge 
excitations in the long-wavelength limit. Generally, one writes the exact response function 
as 

X{k,Uj) = - j— — rr -, (17) 



m(l-G(k,u)) X o(k,uj) ' 
where Xo{k, is the response function 23f] of the noninteracting ID EL and the G(k, u) is the 
charge-charge dynamic local-field factor. In this work, we approximate the dynamic local- 
field factor by its uj— independent, G(k). The charge excitation is determined by calculating 
the pole of response function. The dispersion relation of charge excitations is given by 



u(k) 



A (31-/3* 



1/2 



(A-l) 

where A = exp{nh 2 k / 'mv{k)(l — G{k)) and /3± = h(k 2 /2 ± /c^/c)/m. 



(18) 
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In the limit k — > 0, the effective potential in momentum space, v(k)(l — G(k)) reduce to 
Fourier transformation of bare potential and u(k) turns out to be 



hw(k) = — k J 1 —^ 1 . (19) 



a B 



This behavior of the u(k) is independent of the statistics due to the exchange potential, 
vp(r) is of shorter range than t>s(r) and essentially the low-momentum limit of u(k) remains 
unaltered. The expression can also be written as: 

On the other hand, the charge oscillation eigenmodes provided by the bosonization model 



hamiltonian of the Coulomb interaction liquid 2l|, |24| is given by 



u+ (k) = — ^ + Hi) (1 - Ji + 2mv(k)/nk F ) , (20) 

where g\ = g\kpj1m-K and g\ is a nonsingular interaction matrix element at 2k f- For small 
k excitations, we have 

Q>B y r s 

which is very similar to our finding in Eq. (|T9l) . 



2.3. Fluctuation-dissipation theorem 

In order to calculate the charge excitation spectrum from Eq. (ITS]) , the local-field factor 
is needed. We construct the local-field factor from the fluctuation-dissipation theorem. This 
theorem, which is of paramount importance for systems in equilibrium, relates the dynamic 
susceptibility defined above to the static structure factor 

1 f°° 

S(k) = — / du>$S\x(k,u)) , (22) 
^ Jo 

As x{k, k>) depend on G(k) the above integral expression allows one to determine the local- 
field factor once the static structure factor is calculated. The same approach of obtaining 
the local-field factor has previously been employed by Iwamoto [25] and Dharma-wardana 
and Perrot [26]. We note that the use of fluctuation-dissipation theorem to extract static 
local-field factor is approximate in nature as it neglects the frequency dependence of G(k, u>) 
from the outset. However, the apparent success of our previous work in 2D electron gas 
context |27j also encourages us to use it in the present work. 
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3. NUMERICAL RESULTS 



We turn to a presentation of our numerical results, which are obtained by solving Eq. ([I]) 
with the following self-consistency cycle. We start with the trial choice g(z) = g UF (z) and 
vb(z) = 0, and find the effective potentials v p h{z) by means of Eq. (jSJ) and hence the structure 
factor S(k) via Eq. HI At this point we can calculate new values for g(z) and for vb{z) by 
taking Fourier transforms and using Eqs. (j2J) and (jHJ). This procedure is repeated until self- 
consistency is achieved. We have calculated in this way the pair distribution function of a 
ID EL for various r s and b values. The main results of our work are shown in Figures 1-6. 



Our numerical resu 
DMC data of Ref. 



ts of the pair distribution function of the ID EL is compared with the 
111 ] in Figure 1. The quantum wire width is b = 0.1a* B and the system 
is highly correlated. We have clearly achieved good qualitatively agreement with the DMC 
data. It is obvious from the Fig. 1, FHNC/4+triplet results are in excellent agreement with 
DMC data up to zkp < 1.3 and less agreement in larger z-values to produce the oscillation 
behavior. The strong oscillation behavior in DMC data corresponds to the periodicity of the 
quasi- Wigner crystal. It is apparent from the figure that FHNC/4+triplet approach which 
includes further corrections than FHNC/0, modifies the results and show the first peak 
and its oscillatory behavior, however the magnitude of peak is less than the one predicted 
by DMC simulation. Our g(z) has more structure within FHNC/4+triplet than the g(z) 
obtained within FHNC/0. Moreover, in both approaches the short-range behavior of g(z) 
give correct shape in base on Kimball's cusp condition and they have different values at 
contact. 

In Figure 2, we show our results for S(k) in the paramagnetic ID EL at r s — 1 and 2 within 
both FHNC/0 and FHNC/4+triplet approximations. It is obvious to see the improvements 
brought above by the use of FHNC/4+triplet over FHNC/0. Apart from the strong peak 
structure at q — Akp which appears in DMC simulation, the FHNC/4+triplet results are in 
good agreement with the DMC data of Casula et al [ll| . The strong peak has been related 
to a quasi-order of the electrons [ll|, [2l|. Note that there is no true long-range order in 
ID system. We find that as the density is reduced the correlation effects become stronger 
and S(k) starts to develop a broad peak around q = 4.5&p. The FHNC/0 does not show 
any peak as the density decreases. Theory gives the correct behavior of the long-wavelength 
limit of S(k) and the results coincide with the DMC data. 

There are several noteworthy points based on the results shown in Figs. 1 and 2. (i) 
Despite that FHNC gives very good results for the pair distribution function and static 
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structure factor in high dimensional electron liquids |7|, |8|, [12J, |20j, it can give qualitative 
results for ID EL where the system is highly correlated, (ii) The approach could not produce 
a strong peak for static structure factor at q = 4kp which is related to the slow decay of the 
4kp components of the charge density- density correlation function. 

In order to emphasize the effect of correlations at strong coupling regime, we display 
in Fig. 3 the static structure factor obtained within both FHNC/0 and FHNC/4+triplet 
approximations. The strong coupling arises with decreasing either the density of electrons or 
the quantum wire width parameter b. S(k) move towards right in intermediate of momenta 
value and starts to develop a broad peak around q = 4.5k p within the FHNC/4+triplet 
approach. In the inset of the Fig. 3a, the behavior of the long-wavelength limit of S(k) both 
in FHNC/0 and FHNC/4+triplet approximation are shown in comparison with sum- rule 
given by Eq. (fT41) . As the sum-rule of the long- wavelength limit of S(k) is proven in the 
based on our approach, we numerical evidence that this sum-rule is fulfilled. It is apparent 
from Fig. 3b that by reducing the b value, the peak moves towards right and the magnitude 
of the peak becomes larger for the same r s value. 

The corresponding pair distribution functions are shown in Fig. 4. Again, we show the 
effect of correlations which are corrected by including further correction functions to the 
Bose potential within the FHNC/4+triplet approximation. The value of g(0) decreases with 
decreasing either the density n or b values, but its value always remains positive, consistent 
with the probabilistic definition. Moreover, In the inset of the Fig. 4a, the behavior of 
the short range limit of g(z) both in our approaches are shown in comparison with cusp 
condition. Note that the g(0) value at contact is dependence the approach. We numerical 
evidence that the cusp condition is fulfilled. 

In Fig. 5 we report the effective electron-electron interaction v(z) + vb{z) + vp(z) as it 
emerges from our self-consistent calculations on a ID EL. The attractive part of the effective 
interaction deepens with increased r s or decreased b, as is physically expected. Moreover, 
its deepens in more larger in FHNC/4+triplet than FHNC/0. 

After reaching good agreement between the static structure factor with data from dif- 
fusion Monte Carlo, we can then calculate the correlation energy per particle e c (r s ), the 
difference between the exact ground-state energy and the sum of the noninteracting kinetic 
energy and the exchange energy, of the ID EL as a function of r s , using an integration over 
the coupling strength according to the expression 




(23) 
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(in Rydberg units). Here the static structure factor S\(k) calculated with interaction = 
A(e 2 /e) exp[b 2 k 2 } Ex[b 2 k 2 ]. We have calculated the correlation energy for b = a* B and 2a* B 
over the range 0.1 < r s < 4. The results are reported in Table I in comparison with DMC 
data from Casula et al. [ll|] The Table also include other theoretical results obtained in the 
self-consistent dielectric theory of Singwi et al. 28] with an analytical form of the local-field 
factor from Calmels and Gold. [13( 

It is seen from Table that the present theoretical approach yields fairly better accurate 
values of the correlation energy respect to the Calmels and Gold results [131 ] . 

Finally, in Fig. 6 we display the charge excitation spectrum obtained within the 



FHNC/4+triplet approximation. The u(k) behaviors like k^J\lnk\ at long-wavelength limit 
and reaches to the particle- hole continuum boundary which is given by (3 + in the large fc-limit. 
We use the local-field factor which is calculated from the fluctuation-dissipation theorem by 
Eq. ([22]) . The local- field factor should be useful to study many physical characteristics using 
a suitably defined the dielectric function. 



4. SUMMARY 

In summary, we have presented in this work a theoretical study of the pair distribution 
function, the effective electron-electron interaction and the correlation energy of the ID elec- 
tron liquid where the system is strongly correlated. Our approach yields numerical results 
of good agreement in comparison with recent diffusion Monte Carlo studies|ll|]. We showed 
that the long-wavelength behavior of the static structure factor is in agreement with the 
bosonization findings for the Coulomb Luttinger liquid. The charge excitation spectrum is 
calculated and it is expected to give better results when compared to experimental measure- 
ments in intermediate r s values and even large momentum limit where the Random Phase 
approximation calculation is no longer accurate. Moreover, the small k behavior of the 
charge excitation spectrum is shown be equivalent with the bosonization findings. Improve- 
ments of the theory will be necessary for a quantitative study of the physical quantities 
and the correlation energy. As we have mentioned in the main text, we only considered 
the low-order elementary diagrams and the three-body corrections in the theory is based 
on Fermi hypernetted-chain approximation. We believe inclusion of the contributions of 
the fifth-order and higher elementary diagrams will improve our results. Moreover as we 
mentioned in the main text, we used a practical recipe of FHNC approach [3] by replacing 
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Jastrow product factor to the square of the Slater determinant describing the non-interacting 
fermions gas wave-functions. This recipe of FHNC reduces correlation respect to the full 
FHNC-EL theory. Although it is little affect in high dimensions but it may be important in 
the specific case of a ID EL. 
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FIG. 1: (Color online) The pair distribution function g(z) as a function of zkF for various r s = 1, ( 
top panel) and r s = 2 (bottom panel) at b = 0.1a^_ comparing the FHNC/0 and FHNC/4+triplet 



approximations with DMC data of Casula et al 
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FIG. 2: (Color online) The static structure function S(k) as a function of k/kF for r s = 1 ( top 
panel), r s = 2 ( bottom panel) and 6 = O.la^ comparing both FHNC/0 and FHNC/4+triplet 



approximations with DMC data of Casula et al 
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FIG. 3: (Color online) Top: The static structure function S(k) as a function of k/kp for r s = 3 at 
b = a* B calculating within both FHNC/0 and FHNC/4+triplet approximations. In the inset the 
long- wavelength behavior of S(k) is shown in comparison with sum rule relation (dotted curve). 
Bottom: The static structure function S(k) as a function of kjkp for different r s and b values 
calculating within FHNC/4+triplet approximation. The curves at r s = 2 and 3 have been shifted 
upwards by 0.5 and 1.0, respectively. 
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FIG. 4: (Color online) Top: The pair distribution function g(z) as a function of zfcir for various r s 
values at b = a* B calculating within both FHNC/0 and FHNC/4+triplet approximations. In the 
inset the cusp condition is shown (dotted curve). Bottom: The pair distribution function g{z) as 
a function of zUf for different r s and b values calculating within FHNC/4+triplet approximation. 
The curves at r s = 2 and 3 have been shifted upwards by 0.4 and 0.8, respectively. 
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FIG. 5: (Color online) Top: The total potential v(z) +vb(z) + vf(z) as a function of z/cf for r s = 2 
and 3 values at b = a* B calculating within both FHNC/0 and FHNC/4+triplet approximations. 
Note that deeper potential refer to large r s value. Bottom: The total potential v{z) +vb(z) + vp(z) 
as a function of zkF for various r s at b = a* B calculating within FHNC/4+triplet approximation. 
Deeper potential indicate larger r s value. 
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FIG. 6: (Color online) The charge excitation spectrume u(k) as a function of k/kp for various r s 
values at b = a* B obtaining from Eq. (|22|) within the FHNC/4+triplet approximation. Dashed- 
dotted curve indicates the boundary of particle-hole continuum, f3 + . 
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TABLE I: Correlation energy of the ID EL in Ryd/electron. DMC from Casula, [ll|] STLS from 
Calmels and Gold [l]] 



r s Various calculations b = a* B 



b = 2a* B 



0.1 DMC 

Present work 
STLS 



0.2 DMC 

Present work 
STLS 



0.4 DMC 

Present work 
STLS 



-0.000463 -0.000110 
-0.000459 -0.000107 
-0.000457 -0.000117 



-0.0016996 -0.000418 
-0.001665 -0.000411 
-0.001645 -0.000431 



-0.00579 -0.001514 
-0.00564 -0.001502 
-0.005449 -0.001492 



0.6 DMC -0.01122 -0.003089 

Present work -0.01099 -0.002983 

STLS -0.01044 -0.002955 

0.8 DMC -0.01738 -0.00498 

Present work -0.01678 -0.00476 

STLS -0.01608 -0.00469 

1.0 DMC -0.02394 -0.00709 

Present work -0.02296 -0.00687 

STLS -0.02202 -0.00662 

2.0 DMC -0.05840 -0.01912 

Present work -0.05311 -0.01806 

STLS -0.04968 -0.01735 

3.0 DMC -0.0856 -0.0322 

Present work -0.07330 -0.02952 

STLS -0.06862 -0.02744 

4.0 DMC -0.09986 -0.04372 

Present work -0.08518 -0.03814 

STLS -0.07978 -0.03556 
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